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Abstract 

The integration of the Einstein equations split into the solution of 
constraints on an initial space like 3 - manifold, an essentially elliptic 
system, and a system which will describe the dynamical evolution, 
modulo a choice of gauge. We prove in this paper that the simplest 
gauge choice leads to a system which is causal, but hyperbolic non 
strict in the sense of Leray - Ohya. We review some strictly hyperbolic 
systems obtained recently. 



1 Introduction. 

The Einstein equations equate the Ricci tensor of a pseudo riemannian 4- 
manifold (V, g),of lorentzian signature, the spacetime, with a phenomelogical 
tensor which describes the sources which we take here to be zero (vacuum 
case). The Einstein equations are a geometric system, invariant by diffeo- 
morphisms of V and the associated isometries of g. From the analyst point of 
view they constitute a system of second order quasilinear partial differential 
equations which is over determined, Cauchy data must satisfy constraints, 
and underdetermined, the characteristic determinant is identical to zero. An 
important problem for the study of solutions, their physical interpretation 
and numerical computation is the choice of a gauge, i.e. a priori hypothesis 
for instance on coordinates choice, such that the evolution of initial data 
satisfying the constraints is well posed. 
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The geometric initial data are a 3-manifold M endowed with a riemaniann 
metric and a symmetric 2-tensor which will be the extrinsic curvature of M 
embedded in V = M x R. A natural gauge choice seemed to be the data on 
V of the time lines by their projection N on R and f3 on M. Such a choice 
has been extensively used in numerical computation, though the evolution 
system R^ = obtained with such a choice was not known to be well posed. 
In this article we will show that the system is indeed hyperbolic in the sense 
of Leray-Ohya, in the Gevrey class 7 = 2, and is causal, i.e. the domain 
of dependence of its solutions is determined by the light cone. When the 
considered evolution system is satisfied the constraints are preserved through 
a symmetric first order evolution system. Consideration of a system of order 
4 obtained previously by combination of the equations Ry = with the 
constraints give the same result. 

In section 8 we recall how the old harmonic gauge, interpreted now as 
conditions on N and /3, gives a strictly hyperbolic evolution system and 
the larger functional spaces where local existence and global geometrical 
uniqueness of solutions are known. 

In recent years, since the paper of C-B and York 1995, there has been 
a great interest in formulating the evolution part of Eintein equations as a 
first order symmetric hyperbolic system for geometrically defined unknowns. 
Several such systems have been devised. Particularly interesting are those 
constructed with the Weyl tensor (H. Friedrich, see review article 1996) or 
the Riemann tensor (Anderson, C-B and York 1997) because they lead to 
estimates of the geometrically defined Bel-Robinson energy used in some 
global existence proofs (Christodoulou and Klainerman 1989). We recall in 
the last section the symmetric hyperbolic Einstein-Bianchi system and the 
corresponding Bel-Robinson energy. 

2 Einstein equations. 

The spacetime of general relativity is a pseudo riemannian manifold {V, g),of 
lorentzian signature (- + + +). The Einstein equations link its Ricci tensor 
with a phenomelogical stress energy tensor which describes the sources. They 
read 

Ricci(g) = p 

that is, in local coordinates x x , A = 0, 1, 2, 3, where g = g\^dx x dx^, 
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T3 - pA _ pA , pA pM pA p/i 

where the T's are the Christoffel symbols: 

^ 1 ^ . (9 9 9 . 

The source p is a symmetric 2-tensor given in terms of the stress energy 
tensor T by 

Pap = T a p-^g aP trT, with trT = g x »T XfM 

Due to the Bianchi identities the left hand side of the Einstein equations 
satisfies the identities, with V a the covariant derivative in the metric g 

V Q (-R a/ 3 — -g a p R) = 0, R = g^Rx^ 

The stress energy tensor of the sources satisfies the conservation laws which 
make the equations compatible 

V Q T Q/3 = 

In vacuum the stress energy tensor is identically zero. We will consider 
here only this case. The presence of sources brings up new problems specific 
to various types of sources. 

The Einstein equations (in vacuum) are a geometric system, invariant by 
diffeomorphisms of V and the associated isometries of g. From the analyst 
point of view they constitute a system of second order quasilinear partial 
differential equations which is both undetermined (the characteristic deter- 
minant is identical to zero) and overdetermined (one cannot give arbitrarily 
Cauchy data). 

3 Intrinsic Cauchy problem 

Due to the geometric nature of Einstein's equations it is appropriate to con- 
sider a Cauchy problem also in geometric form. The definition follows. 
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An initial data set is a triple (M, g , K ) where M is a 3 dimensional 
manifold, Qq a riemannian metric on M and K a symmetric 2 tensor. 

An extension of an initial data set is a spacetime (V, g) such that there 

exists an immersion i : M — > Mq C V with i go an d i*Kq equal respectively 
to the metric induced by g on M and the extrinsic curvature of M as 
submanifold of (V,g). 

We say that (V, g) is an einsteinian extension if g satisfies the Einstein 
(vacuum) equations on V . 

A spacetime (V, g) is said to be globally hyperbolic if the set of timelike 
curves, between two arbitrary points is relatively compact in the Frechet 
topology of curves on V. This definition given by Leray 1952 has been shown 
by Geroch 1970 to be equivalent to the fact that (V, g) possesses a Cauchy 
surface, i.e. a spacelike submanifold Mq such that each inextendible timelike 
or null curve cuts Mq exactly once. 

A development is a globally hyperbolic einsteinian extension. 



4 3+1 splitting 

To link geometry with analysis one performs a 3+1 splitting of the Einstein 
equations. We consider a manifold V of the type M x R (the support of a 
development will always be of this type). We denote by x l , G % = 1, 2, 3 local 
coordinates in M, we set x° = t e R. We choose a moving frame on V such 
that at a point (x, t) its axes coincide with the axis of the natural frame, 
tangent to M t = M x {t}, and its axis e is orthogonal to M t , with associated 
coframe such that 9° = dt. A generic lorentzian metric on V with the M t 's 
spacelike reads in the associated coframe 

g = -N 2 dt 2 + gijFP, with 6 i = dx l + (?dt 

The coefficients are time dependent geometric objects on M. The scalar 
N is called lapse, the space vector f3 is called shift, g is a riemannian metric. 
These elements are linked with the metric coefficients g a p in the natural 
frame by the relations: 

.</„ - .</,, n* = (-<p)-\ /r^jvy 

We denote by V the covariant derivative in the metric g. We have 

d d 
di = — , 9 = 9 t - (3 l d h with <9 t = — 
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and we denote 



with Lp the Lie derivative with respect to f3, an operator wich maps a time 
dependent tensor field on M into another such tensor field. 

We denote by K the extrinsic curvature of M t = M x {t} as submanifold 
of (V,g), i.e. we set: 

K ij = ~2N^o9ij 

A straightforward calculation (C-B 1956) gives the fundamental identi- 
ties, written in the coframe 9° = dt, 9 % = dx l + f3 % dt, for the Ricci tensor of 
9 ■ 

Kij = K i:j — ZKjhKi + KijK h — — 

R 0l = N(-V h k? + Vik h h ) 

R m = N{d Kl - NKijK ij + V%N) 



5 Constraints and evolution. 

Constraints 

The following part of the Einstein equations do not contain second deriva- 
tives of g neither first derivatives of K transversal to the spacelike manifolds 
M t . They are the constraints. They read, with S a p = R a /3 — \g a pR, 

Momentum constraint 

Q = ^Roi = -V h K? + ViK* = 

Hamiltonian constraint 

Co = ^5 00 = R- K)K{ + {K h h f = 

These constraints are transformed into a system of elliptic equations on 
each submanifold M t , in particular on M for g= g , K = K , by the con- 
formal method . 
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Evolution. 

The equations 

Kij = Hij — 2K jh K i + K ij K h — — - U 

together with the definition 

dogij = -2NK i:j 

determine the derivatives transversal to M t of g and K when these tensors 
are known on M t as well as the lapse N and shift (3. It is natural to look at 
these equations as evolution equations determining g and K, while N and 
/3, projections of the tangent to the time line respectively on the normal and 
the tangent space to M t , are considered as gauge variables. This point of view 
is conforted by the following theorem (Anderson and York 1997, previously 
given for sources in C-B and Noutchegueme 1988) 

Theorem 1 When R^ = the constraints satisfy a linear homogeneous first 
order symmetric hyperbolic system, they are satisfied if satisfied initially. 

Proof. When = we have, in the privileged frame, 

R = —N 2 R 00 

hence 

S 00 = -R 00 and R = -2N 2 S 00 = 2S° 
2 

and 

S* = AtfiR = -gVS° 

the Bianchi identities give therefore a linear homogeneous system for Sq and 
with principal parts 

N- 2 d S z + g ij djS°, and d S° + d t S l 

This system is symetrizable hyperbolic, it has a unique solution, zero if the 
initial values are zero. The characteristics, which determine the domain of 
dependence, are the light cone. 
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6 Hyperbolicity non strict of = 

An evolution part of Einstein equations should exhibit causal propagation, 
i.e. with domain of dependence determined by the light cone of the spacetime 
metric. 

The equations = are, when N and f3 are known, a second order 
differential system for g^-. The hyperbolicity of a quasilinear system is de- 
fined through the linear differential operator obtained by replacing in the 
coefficients the unknown by given values. In our case for given N, f3 and g^ 
the principal part of this operator acting on a symmetric 2-tensor 7^ is 

\(N~ 2 d 2 00 - g hk 8P hk ) li3 + d% lik + d k d iljk - g h %d jlhk 

The characteristic matrix at a point of spacetime is the linear operator ob- 
tained by replacing the derivation d by a covariant vector £. The charac- 
teristic determinant is the determinant of this linear operator. We take as 
independent unknown 7 12 , 7 23 , 7 31 7 n , 7 2 2> 733 an d consider the 6 equations 
Rij = ,same indices. 

To simplify the writing we compute this matrix in a coframe orthonor- 
mal for the given spacetime metric (N,/3,gij). We denote by (t,x,y,z) the 
components of £ in such a coframe. The characteristic matrix M. reads then 
(up to multiplication by 2): 



'■-z 2 


zx 


zy 





-xy 


xz 


t 2 -x 2 


xy —yz 








yz 


xy 


t 2 -y 2 


—xz 





2xy 





2zx t 2 -y 2 - z 2 


-x 2 


-x 2 


2xy 


2yz 


-y 2 


t 2 -x 2 - z 2 


-y 2 





2yz 


2xz —z 2 


-z 2 


t 2 — x 2 — y 2 



The characteristic polynomial is the determinant of this matrix. It is 
found to be 

DetM = b 6 a 3 , with b = t, a = t 2 - x 2 - y 2 - z 2 

The characteristic cone is the dual of the cone defined in the cotangent plane 
by annulation of the characteristic polynomial. For our system the character- 
istic cone splits into the light cone of the given spacetime metric and the 
normal to its space slice. Since these charateristics appear as multiple and 
the system is non diagonal it is not hyperbolic in the usual sense. We will 
prove the following theorem 
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Theorem 2 When N > and (3 are given, arbitrary, the system Rij = 
is a system hyperbolic non strict in the sense of Leray Ohya for g^, in the 
Gevrey class 7 = 2, as long as gij is properly riemannian. If the Cauchy data 
as well as N and f3 are in such a Gevrey class the Cauchy problem has a local 
in time solution, with domain of dependence determined by the light cone. 

Proof. 

The product ab 2 M.~ l , with .M -1 the inverse of the characteristic matrix 



M. is computed to be: 










t 2 — x 2 — y 2 —zx 


-zy 








xy 


—zx t 2 — y 2 — z 2 


-xy 


zy 








-zy -xy 


t 2 -x 2 - z 2 





zx 





-2xy 


—2zx 


t 2 -x 2 


x 2 


x 2 


—2xy —2zy 





y 2 


t 2 -y 2 


y 2 


-2zy 


—2zx 


z 2 


z 2 


t 2 -z 2 



We see that the elements of the matrix ab 2 M.~ x are polynomials in x, y, z. 
The product of this matrix by Ai is a diagonal matrix with elements ab 2 in the 
diagonal. Consider now the differential operator acting on g^. Multiply 
it on the left by the differential operator defined by replacing in ab 2 M.~ x the 
variables x,y,z by the derivatives 81,82,03. The resulting operator is quasi 
diagonal with principal operator 8 X 8\8 2 . It is the product of two strictly 
hyperbolic operators, 8 x 8\8q and 8q. The result follows from the Leray-Ohya 
general theory. 

7 Hyperbolic non strict 4th order system. 

Lemma 3 The following combination of derivatives of components of the 
Ricci tensor of an arbitrary spacetime : 

Ay = dodoRij — doV(iRj)o + Vj<9ji?oo 

reads, when g is known, as a third order quasi diagonal hyperbolic system for 
the extrinsic curvature K~ij. 

Aij = d VK ij + d d (HK ij -2K im K™)-d d (N- 1 V j d i N) + 

d (-V (i (K j)h d h N) - 2NR h :ijm K™ - NR m(t K% + HVfrN) + 
Vjdi(NAN - N 2 K.K) + 
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with 



VK .. = -80^-%^) + V h V h (NK tJ ), A = V h V h , H = K* 

and 

dj = VjdiiNdoH) - d (NVjdiH) 

Proof. A straightforward computation shows that Cij contains terms of 
at most second order in K (and also in N) and first order in g (replace d^g^ 
by —2NKij). The other terms of A^-, except for d VKij are second order 
in K. All terms of A^- are at most second order in g except for third order 
terms appearing through Vjdi(NAN). Because of these terms the system 
for g and K given by 

A tJ = (1) 

and 

d o9ij = -2NKij (2) 

is not quasi diagonal. It is not hyperbolic in the usual sense of Leray, we will 
prove the following theorem. 

Theorem 4 The system (1), (2) with unknown g, K is for any choice of 
lapse N and shift f3 equivalent to a system hyperbolic non strict in the sense 
of Leray-Ohya with local existence of solutions in Gevrey classes 7 = 2 and 
domain of dependence determined by the light cone. 

Proof. Replace in the equations A^ = the tensor K by — {2N)~ l dog: 
this gives a quasi diagonal system for g, but with principal operator (do) 2 d x d\. 
The result follows immediately from the Leray-Ohya theory. 

The system for g, K can be turned into a hyperbolic system by a gauge 
choice as follows. 

Theorem 5 Suppose that N satisfies the wave equation. 

N- 2 d d N — AN — f (3) 

with f an arbitrarily given smooth function on MxR. The system (1),(2),(3), 
called S, is equivalent to a hyperbolic Leray system for g, K, N, for arbitary 
shift. 
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Proof. We use the wave equation (3) to reduce the terms — dQ (N~ l V jdiN) + 
Vjdi(NAN) in A^- to terms of third order in N, second order in g and K. 
We replace the equation (3) by the also hyperbolic equation 

d (N~ 2 d d N- AN) =d f (4) 

and replace dog by —2NK wherever it appears. Then the equation (4) is 
third order in N, while first order in g and K. We call S' the system thus 
modified. 

We give to the equations and unknowns the Leray-Volevic indices: 

m(l) = 0, m(2) = 2, m(3) = 1 (5) 

n(g) = n(K) = 3, n(N) = 4 (6) 

The principal matrix of the system S' is then diagonal with elements the 
hyperbolic operators dod a d a or do. 

If the equation (3) is satisfied on the initial submanifold as well as S', the 
equation (3) and the system S are satisfied. 

Remark 6 The system A^ = has the additionnal property to satisfy a 
polarized null condition, that is the quadratic form defined by the second 
derivative of A at some given metric g, A" ij(g)(j, 7), vanishes when 7 = t®t 
with £ a null vector for the spacetime metric g such that 7 is in the kernel of 
the first derivative of the Ricci tensor of spacetime at g ( C-B 2000). 

8 An hyperbolic second order system 

A variety of hyperbolic evolution systems for Einstein equations have been 
obtained, with a speed greatly increasing in recent years, by replacing the 
trivial gauge choice (which is the data of and (3 on V) by more elaborate 
ones, together with combining the evolution equations with the constraints. 
The hope in changing the gauge is to find systems either better suited to the 
study of global existence problems, or more stable under numerical codes. We 
give some references in the bibliography. We will return below to the original 
gauge choice (C-B 1952) in the perspective of conditions on the lapse and 
the shift. 

The following identity was already known by De Donder, Lanczos and 
Darmois. It splits the Ricci tensor with components R a p in the natural 
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frame into a quasi linear quasi diagonal wave operator and 'gauge' terms, as 
follows: 

with 



where H a/3 is a quadratic form in first derivatives of g with coefficients poly- 
nomials in g and its contravariant associate. 
The F x are given by 

F x = = V a V a x w 

They are non tensorial quantities, result of the action of the wave operator 
of g on the coordinate functions. For this reason their vanishing is called 
'harmonicity condition'. 

The contravariant components of the Ricci tensor admit an analogous 
splitting, namely: 

pa/3 ^ r>af3 1 / a \ 9 n o X 8 

with 

1 d 2 

R *P = _ A M a/3 + R aP 

Let us denote by R K a » the components of the Ricci tensor in the previous 
frame 6° = dt, 9 l = dx % + (3 l dt ) to distinguish them from the components in 
the natural frame now denoted Rap. These components are linked by the 
relations: 

w g(dg) djdx^) _ 



doo _ pap djdt) d(dt) _ 00 

- H w W W ~ {e) 

P 0i _ pa/3 9(6 ) d{6 % ) _ oo R oj 

{e) d{dx a ) d{dxP) 

The equations = are a quasidiagonal second order system for g^ when 
N and (5 are known. The equations R^ = and are quasilinear wave 
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equations for N and (3 when the g'^s are known: we interpret these equations 
as gauge conditions. The set of all these equations constitute a quasidiagonal 
second order system for g^, N and (3, hyperbolic and causal as long as N > 
and g is properly riemannian. 

The Bianchi identities show that for a solution of these equations the 
quantities F x satisfy a linear homogeneous quasidiagonal hyperbolic and 
causal system. Its initial data can be made zero by choice of initial coor- 
dinates if and only if the geometric initial data g, K satisfy the constraints. 
A solution of our hyperbolic system satisfies then the full Einstein equations. 

The following local existence and uniqueness theorem improves the dif- 
ferentiability obtained in the original theorem of C-B 1952 who used C k 
spaces and a constructive (parametrix) method. The improvement to Sobolev 
spaces (one can endow M with a given smooth riemannian metric to define 
those spaces) with s > 4 for existence and s > 5 for geometric unique- 
ness is given in C-B 1968 using Leray's results. The improvement given in 
the theorem was sugested by Hawking and Ellis 1973, proved by semigroup 
methods by Hughes, Kato and Marsden 1978 and by energy methods by 
C-B, Christodoulou and Francaviglia 1979. The other hyperbolic systems 
constructed in the past twenty years did not lead, up to now, to further 
improvement on the regularity required of the Cauchy data. Such an im- 
provement would be an important step. 

Theorem 7 Given an initial data set, g~o, Kq G H l s ocal , H l s °ff satisfying the 
constraints, there exists an einsteinian extension if s > 3. 

The question of uniqueness is a geometrical problem. It is in general 
easy to prove that the solution is unique in the chosen gauge, for instance in 
the harmonic gauge recalled above. But two isometric spacetimes must be 
considered as identical. The following theorem (C-B and Geroch) gives this 
geometric uniqueness (maximal means inextendible). 

Theorem 8 The development of an initial data set is unique up to isome- 
tries in the class of maximal developments if s > 4. The domain of depen- 
dence is determined by the light cone of the spacetime metric. 

The proof of C-B and Geroch considers smooth data and developments, 
the refined result is due to Chrusciel 1996. The geometric uniqueness in the 
case s = 3, even the local one, is still an open problem. 
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9 Bianchi equations. 

The Riemann tensor satisfies the identities 

(7) 
it holds therefore that, modulo the symmetries of the Riemann tensor 

^<xR\.H,Py = ^0R~ft* ~ (8) 

hence if the Ricci tensor R a/3 satisfies the vacuum Einstein equations 

RaP = (9) 

it holds that 

V Q it^ 7 = 0. (10) 

The system (8), (11) splits as the Eintein equations into constraints, con- 
taining no time derivatives of curvature, namely in the frame used in the 3+1 
splitting: 

^iRjk,\n + VkRijM + V ' jRki,\n = (11) 
V a R\^ = 0. (12) 

and an evolution system 

Vo-Rm;,a,u + Vfc-R /i,A/i — Vfe-R fc,A^ = (13) 
+ V h Rt, M = (14) 
This system has a principal matrix consisting of 6 identical 6 by 6 blocks 
around the diagonal, obtained by fixing a pair A,/i, A < //. Each block is 
symmetrizable through the metric g, and hyperbolic if g is properly Rieman- 
nian and N > because the principal matrix M° for the derivatives do was, 
up to product by N^ 1 , the unit matrix and the derivatives dh do not contain 
d/dt. 

The Bel-Robinson energy is the energy associated to this symmetric hy- 
perbolic system. 

Remark 9 Following Bel one can introduce two pairs of gravitational "elec- 
tric" and "magnetic" space tensors associated with the 3+1 splitting of the 
spacetime and the double two-form Riemann(g) : 

N 2 E t] ee Ro h0v D.^-^^R^ 

1 hk 1 hk 

NHij = -^rj ihk R..... ^ NBji = -r] ihk A$^ 
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where r] i j k is the volume form of g. The principal part of the evolution sys- 
tem ressemble then to the Maxwell equations, but contains an additional non 
principal part. Its explicit expression is given in Anderson, C-B and York 
1997. 

The Bianchi equations do not tell the whole story since they contain the 
spacetime metric g, which itself depends on the Riemann tensor. 

A possibility to obtain a symmetric evolution system (Friedrich 1996 with 
the Weyl tensor) for both g and Riemann(g) (Anderson, C-B and York 1997) 
is to introduce again the auxiliary unknown K and use 3 + 1 identities, in- 
volving now not only the Ricci tensor but also the Riemann tensor. One can 
then obtain a symmetric first order hyperbolic system for K and T, the space 
metric connection, modulo a choice of gauge, namely the integrated form of 
the harmonic time-slicing condition used before (C-B and Ruggeri 1983). 
The energy associated to this system has unfortunately no clear geometri- 
cal meaning. Determination of the metric from the Riemann tensor through 
elliptic equations seems more promising for the solution of global problems 
(see Christodoulou and Klainerman 1989, Andersson and Moncrief, in prepa- 
ration) 
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